A relativistic partially electromagnetic planar plasma shock 
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We model relativistically colliding plasma by PIC simulations in one and two spatial dimensions, 
taking an ion-to-electron mass ratio of 400. Energy dissipation by a wave precursor of mixed polarity 
and different densities of the colliding plasma slabs results in a relativistic forward shock forming 
on millisecond timescales. The forward shock accelerates electrons to ultrarelativistic energies and 
reflects upstream ions, which drag the electrons along to preserve the plasma quasi-neutrality. No 
reverse shock forms. The shock may be representative for internal gamma ray burst shocks. 
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PACS numbers: 98.70.Rz,52.35.Tc,52.65.Rr,52.27.Ny 



Accreting compact objects like neutron stars or black 
holes are amongst the most energetic environments in 
the universe. Such systems can range from micro- 
quasars through active galactic nuclei to gamma ray 
bursts (GRBs) 0, H, IE El- Common to these objects 
is the ejection of a relativistic jet that is probably accel- 
erated by magnetic fields [EJ(|. GRB jets are ejected by 
a forming compact object [7(. Their high Lorentz factor 
100 < r < 10 3 introduces a radiation beaming, making 
GRBs detectable at cosmological distances. The GRB 
emissions are explained either by the fireball model, in 
which the plasma kinetic energy dominates EL Q , or by 
the electromagnetic model, in which the Poynting flux is 
dominant [8j, [9( . We consider the fireball model. 

The high energies involved in the thermalisation of rel- 
ativistic jets imply that plasma structures form, such 
as phase space hole s flCj . A shock description needs 
a kinetic approach 11]. Particle-in-cell (PIC) simula- 
tions have modelled the external GRB shock between 
the jet and the ambient plasma El El El EE EE E3 , 
and have examined the magnetic field amplification and 
particle acceleration. The weakly magnetized external 
GRB shocks develop out of a broad wave spectrum 
El El, HE HI , and are thus filamentary. The inter- 
nal shocks, powering the prompt emissions, form closer 
to the collapsing star with its strong magnetic fields. The 
magnetic field pressure may initially exceed the ther- 
mal pressure of the j et p lasma, in contrast to the case 
considered by Ref. B3- The magnetic field and the 
high plasma temperature [22j suppress the filamentation 
modes [23l . l24l |25| and other wave modes can develop. 

The finite grid instability requires a high spatial resolu- 
tion on electron timescales of multi-dimensional PIC sim- 
ulations [III, while a relativistic shock will extend over 
ion scales. A large ion-to-electron mass ratio is required 
to model shocks [27|- To account for these constraints, 
we select the initial conditions such as to physically con- 
fine the plasma dynamics to one spatial dimension (ID) 
and three momentum components, which allows us to re- 
strain also the simulation geometry to one spatial dimen- 
sion. We verify the validity of this assumption with the 
help of a two-dimensional (2D) simulation. We demon- 



strate that the magnetic field results in the formation 
of a relativistic planar shock. The energy dissipation is 
provided by a wave with mixed electrostatic and electro- 
magnetic polarity, probably a nonlinear oblique whistler 
[28| . This contrasts the energy dissipation by the purely 
electromagnetic filamentation modes considered in Refs. 
EE El EE EE EE SI, and the nonrelativistic electro- 
static shocks modelled in Refs. [IE HH- Different slab 
densities [29| facilitate the development of a shock with 
a Mach number of 40 and an AlfVenic Mach number of 
10. The forward shock reflects upstream ions. 

PIC simulation methods [3(| solve the Maxwell equa- 
tions and the relativistic equation of motion 
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and fulfill V • B = and V • E = p/eo as constraints. The 
electric E and magnetic B fields are defined on a discrete 
grid. Each computational particle (CP) with index i of 
the species j has the charge qj and mass rrij with qj / rrij — 
q/m, where q and m are the physical charge and mass of 
the particles, represented by the species j. The particle 
momentum is Pi = mjVjr(Vj). 

Two plasma slabs collide in the simulations at the po- 
sition x = 0, each consisting of the electrons and ions. 
The ion-to-electron mass ratio is R = rrii/m e = 400. The 
plasma frequency of the species j with the number den- 

1/2 

sity rij is uj Pt j = (rijq'j /mjEo) . The electrons of slab 

1 have = 10 5 x 2n/s and the ions uj p2 = lo P} \/\/R~. 
The dense slab 1 is initially moving with the positive 
speed Vb along the x— direction. The electrons and ions 
of slab 2 have ui p3 = uj p i/\/To and lu p a — u} p ^/yR, re- 
spectively. The tenuous slab 2 moves with the speed —Vf, 
along the x— direction. Initially, all mean speeds along 
the y, z directions vanish. All plasma s pec ies have the 
observed jet temperature T = 100 keV [22]. The posi- 
tion and time have the units of the electron skin depth 
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FIG. 1: (Color online) The electric fields in V/m in the 2D 
simulation at t = ti'- (a) shows the flow-aligned electric field 
component E x and (b) shows E y , the second electric field 
component in the simulation plane. 
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FIG. 2: (Color online) The 10-logarithms of the electron dis- 
tributions in the ID simulation at t = t%: The distribution 
in the x,p x plane is shown in panel (a). Panels (b) and (c) 
display the distributions in the x,p y and x,p z planes. 



A. s 



c /«: 



>,3 



) 1 ^ 2 and Aj 



/c, where c is the speed of 



light in vacuum. Each slab initially occupies half of the 
simulation box. The slabs collide along the x— direction 
at the speed v c — 2Ufe/(l + v$) = 0.9 c, which is rep- 
resentative of internal GRB shocks @. We introduce 
a spatially homogeneous oblique magnetic field. The 
magnetic field component with eB x /m e = s/2u) p> i sup- 
presses the filamentation instability [23| and a perpendic- 
ular B z = B x /10 enhances the energy dissipation. Ini- 
tially, we set By = 0. The initial magnetic field strength 
is thus 5 (J.T. The initial box-averaged plasma kinetic 
energy density exceeds the magnetic field energy den- 
sity by a factor of 10 4 . The ID simulation box length is 
L = 1.58 x 10 4 A S , which is resolved by 1.8 x 10 5 cells. 
The 2D box has the extent L x = 714 X s (1.6 x 10 4 cells) 
and L y = 11 A s (200 cells) along the x and y directions, 
respectively. The boundary conditions are periodic. The 
ID (2D) simulation time is h = 4370 (t 2 = 194). The 2D 
simulation represents each plasma species by 64 particles 
per cell (PPC). The ID simulation resolves each species 
of the dense (tenuous) slab by 576 (256) PPC. 

The electric field in the 2D simulation plane at t = t 2 in 
Fig. [1] reveals a planar structure at x ps 120, the leading 
edge of the dense plasma slab. A similar E z distribution 
is not shown. The filamentation instability would lead to 
the structure formation on a scale of X s along y, and it 
is here suppressed by B x [23|. No strong field is present 
at x ps —120. We exploit the planarity of the slab front 
in Fig. Q] and resort to the ID simulation with its better 
statistical plasma representation. 

The current system responsible for the electric fields 
in Fig. Q] is displayed in Fig. [2] The electrons of the 
tenuous slab are deflected by the oblique magnetic field 
that is frozen into the dense downstream plasma. Since 
the massive upstream ions are not significantly deflected, 
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FIG. 3: (Color online) The 10-logarithms of the electron dis- 
tributions in the ID simulation at t — 1550: Panel (a) shows 
the distribution x,p x . The distributions in x,p y and x,p z are 
exhibited in panels (b) and (c), respectively. 



a current and thus the electric field is building up at 
the leading edge of the plasma slab in Fig. [1] Further 
downstream, Fig. [Ha) shows the electron phase s pac e 
holes, which thermalize the electrons, as in Ref. [25j |. 
These are also observed at the Earth bow shock [3lj . 

The energy dissipation by the electric fields eventually 
modulates the electrons of the dense plasma slab, giving 
rise to larger electric fields. Consequently, Fig. [3] demon- 
strates that the electrons are accelerated to T factors of 
tens to hundreds and follow a corkscrew orbit. The rel- 
ativistic mass of the electrons is not small compared to 
mi = 400 m e at t = 1550. Figure Ufa) shows the ion re- 
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FIG. 4: (Color online) Panel (a) shows the ion distribution in 
x,p x space. The color is the 10-logarithm of the number of 
CPs. Panel (b) shows the electrostatic E x field and panel (c) 
shows E y and cB z (red). The time is t = 1550. 
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FIG. 5: The shock: (a) shows the ion distribution in x,p x , 
while (b) shows the corresponding electron distribution. The 
color is the 10-logarithm of the number of CPs. 



sponse to the strong fields. The spatial interval with the 
strongest particle acceleration at x ss 780 is well behind 
the expected slab front vt, t/X s ss 980 at t = 1550. Some 
ions have been accelerated beyond the initial velocity Vi, 
by the electric field at the slab front [25| and they have 
propagated farther upstream beyond x ss 980. 

The wave modes at x w 800 that grow only close to 
the leading edge of the dense plasma slab, as Fig. BJb) 
exemplifies for the E x field component, could have origi- 
nated from an oblique whistler instability, similar to that 
discussed in Ref. [28|]. The wave growth at x < is sup- 
pressed or delayed. This will introduce an asymmetric 
plasma dynamics in the intervals separated by x = 0. 

The wave amplitudes in the considered plasma scale 
linearly with io v .\ if all ratios between w p j remain the 
same as discussed, for example, in Ref. [32j . The 
wave is circularly polarized, which explains the electron 
corkscrew orbits in Fig. [3l The peak value of B z in 
Fig. |3Jc) exceeds the initial background magnetic field 

B = {B 2 X + B 2 Z ) 1/2 with cB a 1500 V/m, and the wave 
in Fig. |3Jc) is nonlinear. The mixed polarity of the 
waves allows them to interact with the plasma through 
their electrostatic component. The wave amplitude is 
sustained, because the wave is driven by the dense plasma 
slab that compensates damping losses. The damping may 
explain why E y and B z are in phase in Fig. [4jc). 

At t = <i, the ion distribution in Fig. Ufa) has de- 
veloped into a shock. The foreshock in the interval 
2000 < x < 3200 is formed by incoming and shock- 
reflected ions. The electron transport with the shock- 
reflected ions across the upstream B in Fig. 0(b) en- 
sures the plasma quasi-neutrality. The two ion beams 
at x < 500 result from the delayed shock formation. At 
x 1800, the p x values of the electrons are comparable 
to those of the ions multiplied by mi/m e in the simu- 



lation reference frame. A Lorentz transformation of the 
electron velocity into the upstream frame, which moves 
with v = —vt, relative to the simulation box frame, yields 
a peak electron T ~ 700 for our v c = 0.9 c. In com- 
parison, the filamentation instability in Ref. [l7[ yields 
only a peak electron Y w 10 for a collision speed v/ with 
Y(vf) — 5. The planar electromagnetic shock considered 
here is a stronger electron accelerator. 

The shock in Fig. [5] has formed during the time t = ti, 
which is 7 ms in physical units. If we assume that this 
shock would constitute an internal shock of a GRB and 
that the jet frame moves at a Y ss 10 3 , the time dilation 
would yield a duration of seconds for the shock develop- 
ment in the Earth (observer) frame of reference. This 
is comparable to the timescale of the prompt emissions. 
The shock we consider is likely to form during a longer 
time for rrii/m e = 1836. The bulk Lorentz factor of GRB 
jets may, however, also be less than 10 3 @. Different jet 
plasma densities, which may be higher than we consider 
in this work [33|], do not strongly affect the shock devel- 
opment time due to a scaling of the timescales with the 
square root of the density [32[ . 

In summary, we have examined the collision of two 
plasma slabs with different densities at a relativistic 
speed. The bandwidth required to model the plasma 
dynamics at a high ion-to-electron mass ratio cannot be 
realized at present with PIC simulations in more than 
one spatial dimension due to the necessary large box 
size [3j and the high-resolution grids that are required 
in higher dimensions to suppress the electromagnetic fi- 
nite grid instability [261 ] . The quasi-parallel magnetic 
field has, however, suppressed the beam filamentation, 
and the initial planar shock development may thus be 
well-approximated by a one-dimensional PIC simulation. 
The energy dissipation has been accomplished through a 
nonlinear wave of mixed polarity, which could have origi- 



4 



nated from an oblique whistler instability. Whistlers are 
frequently observed also at slower quasi-parallel plasma 
shocks and they are known to be efficient particle ac- 
celerators [H]]. The large energy dissipation, which ulti- 
mately forces the shock to form, is likely to distort the 
guide magnetic field and the slab boundary would likely 
be deformed by the ram pressure excerted by the up- 
stream plasma. However, a significant distortion of the 
guide magnetic field away from the flow-aligned direction 
introduces a stronger perpendicular magnetic field com- 
ponent. The latter will probably not suppress or delay 
the shock formation and the plasma thermalization. Our 
simulation may thus give a correct estimate for the shock 
formation time even in 3D. 

The shock we consider here may thus be representative 
for a GRB internal shock. Such a rapid shock formation 
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